QUANTUM CHAOS IN THE HEAVY QUARKONIA 
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The binding energy spectra of the heavy quarkonia are calculated by solving the Schrodinger 
equations with Coulomb plus confining potentials. Statistical properties of the obtained spectra are 
studied by plotting nearest level spacing distribution histograms. 

PACS numbers: 05.45 Mt; 14.40.Lb; 14.40.Nd; 12.39.Pn 



I. INTRODUCTION 

Level statistics, i.e. the statistical theory describing 
the fluctuation properties of the quantum spectra was 
originally developed in nuclear physics in early fifties. 
Pioneering work dates back to Wigner p| who gave the 
first detailed theoretical analysis of the nuclear spectra 
from the viewpoint of random matrix theory. Further de- 
velopments of this subject were carried out by Porter , 
Mehta Q and Dyson Q. In particular, it has been estab- 
lished that the spectra of the complex quantum many- 
body systems have a behaviour similar to that of the 
random matrices 

Later the development of the statistical spectroscopy 
was stimulated by quantum chaology which arose in early 
cighthies as a quantum theory of classical chaotic sys- 
tems. It was established that quantized spectra of classi- 
cal chaotic system spectra have the same statistical prop- 
erties as those of the random matrices. It is well known, 
that the classical chaotic dynamics is characterised by 
an exponential divergence of the neighboring trajectories 
in phase-space. How far the complexity of the classical 
motion is reflected in the wave dynamics of the corre- 
sponding quantum system (described by a linear wave 
equation) is an interesting aspect of the study of deter- 
ministic non-linear dynamics. Numerical evidence has 
shown that the statistical properties of the energy spectra 
of the quantum systems depend crucially on the degree of 
chaos in the underlying classical dynamics. Namely, the 
spectrum of the quantum operator, describing classically 
non-integrable systems, has the same statistical proper- 
ties as those of random matrices. It has been established 
that certain distribution functions and fluctuation mea- 
sures for a wide variety of complex nuclear spectra fol- 
low a universal behavior and that moreover they agree 
very well with the results obtained by ensemble aver- 
aging over certain matrix ensembles with randomly dis- 
tributed matrix elements. In particular, for systems with 
time-reversal symmetry, good agreement was found with 
a Gaussian orthogonal ensemble(GOE) [(Jut- It has also 
been established, that for integrable systems the nearest 
neighbor spacing distribution is dominated by level clus- 



tering and is expected to be Poisson-like. Later, this has 
been confirmed in a number of studies of model systems 
like billiard balls, simple rotators etc.[(| M, as weu as 
in realistic systems such as hydrogen atom in a magnetic 
field GUim. 

Details of the progress of level statistics and 
random matrix theory are given in a recent review 12] . 
Here we deal with the spectral statistics of the quark- 
anitquark systems (quarkonia). Quarkonium spectrosopy 
has been the subject of extensive experimental as well 
as theoretical studies ^3|"E3- Theoretical description 
of quarkonium properties can be done quantum mechan- 
ically as well as within the QCD approach. Quantum 
mechanical description of hadrons is reduced to solving 
wave equations with Coulomb plus a confining potential. 
Considerable progress has been made in the calculation of 
the spectra of quarkonium an d baryons ^tJ- At 

present a large amount of theoretical and experimental 
data is available. Recently energy fluctuations and quan- 
tum chaos in hadrons and QCD has become a subject of 
theoretical studies [l8|-[24|. In partuclar, it has been 
found that QCD is governed by quantum chaos in both 
confinement and deconfinement phases [l8L I23I ] . The sta- 
tistical analysis of the measured meson and baryon spec- 
tra shows that there is quantum chaos phenomenon in 
these systems |20j . The study of the charmonium spectral 
statistics and its dependence on color screening has es- 
tablished quantum chaotic behaviour of this system [24| . 
It was claimed that such a behaviour could be the reason 
for J supression [24|. 

In the present work we concentrate our attention on 
level fluctuations and quantum chaos in heavy quarko- 
nia due to the confining forces. It is well known that 
many realistic and model confined dynamical systems 
such as an atom in a magnetic field, particle motion in 
resonators, quantum dots, billiards exhibit chaotic dy- 
namics. Quarks in the hadrons are also confined sys- 
tems, and to some extent are equivalent to motion in 
billiards (domain with hard walls where elastic scattering 
of the particle will occur ) [2M l26l Wfy . Depending on the 
shape of the billiard domain the motion can be regular 
or chaotic [2g. If one uses for quark- antiquark poten- 
tial Coulomb plus harmonic oscillator potential, the sys- 
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tern becomes equivalent to the atom in a magnetic filed, 
which is also described by a Coulomb plus harmonic os- 
cillator potential whose dynamics is chaotic in classical 
and quantum cases [jjj, [TlJ • An analysis of the quarko- 
nium would suggest that the dynamics of such a system 
could be chaotic in classical as well as in quantum cases 
. Therefore the present paper may be viewed as an at- 
tempt to use well-known technique to the quarkonium 
systems to establish their chaotic behaviour. 

Charmonium and bottomonium spectral statistics 
based on the solution of the Schrodinger equation for 
this system is considered . Two types of potentials 
for quark-antiquark interaction are considered: Coulomb 
plus linear and Coulomb plus harmonic oscillator poten- 
tials. Solving the Schrodinger equation with these po- 
tentials by the scaling variational method we obtain the 
energy levels of the quarkonium system. This paper is 
organized as follows: In the next section we solve the 
Schrodinger equation analytically using the scaling vari- 
ational method and find the mass sectra of charmonium 
and bottomonium. In section 3 we calculate the statis- 
tical properties of the spectra of charmonium and bot- 
tomonium and in section 4 some conclusions and final 
remarks are presented. 



II. SCHRODINGER EQUATION FOR 
QUARKONIUM 

To obatain the quantum spectra the following 
Schrodinger equation is solved 
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where Z = |a s , with a s being the strong coupling 
constant, A is the parameter of the confining potential, 
k = 1, 2 and M = m q /2 is the reduced mass of quarko- 
nium (for simplicity we consider quarks with the same 
flavor). To solve this equation we use the scaling vari- 
ational method. This method was shown to work 
well for Coulomb plus confining type potentials. Recently 
this method was used for the calculation of doubly heavy 
baryon spectra [2!j and heavy-light meson spectra |30f . 
For analysis, eq. Q may be written as 



where 



and 



(H + H')^ ni = E nl * nl , 



H = -U-Z, 
2 r 



H' = Xr k 



(2) 



(3) 



(4) 



According to the scaling variational method we choose 
the trial wave function of the form 



where 4> n i(br) are the the eigenfunctions of the Hamil- 
tonian Hq, i.e. 

HQ(j) n l — e n i4> n i, 
satisfying the the normalization condition 

< 4>nl\4>n'V >= S nn 'S W . 

and the eigenvalues are given as 



2n 2 



Then the eigenvalues of eq.QJ can be calculated as 

E n l = £nl(b ), (6) 

where 

e n l(b) = b 3 < cj>(b,r)\H + H'\4>{b,r) >, 
and &o is found by minimizing this energy, i.e. 
de 
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0. 



(7) 



The quantity e n i(b) may be calculated as: 



h 2 nh~ k 
e nl (b)=n- 2 C--Zb-^-), 



where 



q = kXn 2 < <p n i\r k \<f>„i > 



These matrix elements are calculated analytically [31f . 
For k = 1 we have 

q = Xn 2 ^(3n 2 + 

where n and / are the principal and orbital quantum num- 
bers, respectively. Then the extremum condition gives 
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* ni (6,r) = b^ 2 (f> nl (br), 



(5) 



3 3(27 g + 2Z 3 + 3%(27g + 4Z 3 ))2) 

3(27q + 2Z 3 + 2,i{q{27q + AZ 3 ))^ 
2^3 

For k = 2 we have 

g = 5n 2 + l -3/(7 + 1). 

The extremum condition JJJ leads to an equation for 
defining b 

bf ) -Zb 3 -q = 0. 

Solving this equation we get bo that leads to the en- 
ergy eigenvalues for the case k — 2. Thus the solutions of 
the Schrodinger equation for Coulomb plus confining po- 
tentials are obtained analytically using the scaling varia- 
tional method. 
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III. SPECTRAL STATISTICS 

Now an analysis of the statistical properties of the 
spectra is to be considered. One of the main charac- 
teristics of the statistical properies of the spectra is the 
level spacing distribution 0,0,13 function. We calcu- 
late the nearest-neighbor level-spacing distribution for 
quarkonium spectra. The nearest neighbor level spacings 
are defined as Si = Ei+i — E%, where E{ are the energies 
of the unfolded levels, which are obtained by the follow- 
ing way: The spectrum {Ei\ is separated into smoothed 
average part and fluctuating parts. Then the number of 
the levels below E is counted and the following staircase 
function is defined: 

N(E)=N av {E)+N fluct {E). 

The unfolded spectrum is finally obtained with the 
mapping 

Ei = N av (Ei). 

Then the nearest level spacing distribution function 
P(S) is defined as the probability of S lying within the 
infinitesimal interval [S, S + dS] . 

For the quantum systems which are chaotic in the clas- 
sical limit this distribution function is the same as that of 
the random matrices [|| 0] . For systems which are reg- 
ular in the classical limit its behaviour is close to a Pois- 
sonian distribution function. This distribution is usually 
taken to be a Gaussian with a parameter d: 

P(H) ~ exp(-Tr{HH+}/2d 2 ), 

and the random matrix ensemble corresponding to this 
disribution is called the Gaussian ensemble. For Hamilto- 
nians invariant under rotational and time-reversal trans- 
formations the corresponding ensemble of matrices is 
called Gaussian orthogonal ensemble (GOE). It was es- 
tablished 0, |U that GOE describes the statistical 
fluctuation properties of a quantum system whose clas- 
sical analog is completely chaotic. The nearest neighbor 
level spacing distribution for GOE is described by Wigner 
distribution: 

P(S) = ^Sexpi—nS 2 ). (8) 

The common way to study of the level statistics is to com- 
pare the calculated nearest-neighbor level-spacing distri- 
bution histogram with the Wigner distribution. 

For systems whose classical motion is neither regular 
nor fully chaotic (mixed dynamics) the level spacing dis- 
tribution will be intermediate between the Poisson and 
GOE limits. Several empirical functional forms for the 
distribution were suggested for this case ^ the 

Hamiltonian is not time-reversal invariant, irrespective 
of its behavior under rotations, the Hamiltonian matrices 
are complex Hermitian and the corresponding ensemble is 



called Gaussian unitary ensemble (GUE). If the Hamilto- 
nian of the systems is time-reversal invariant but not in- 
variant under rotations, then the corresponding ensemble 
is called the Gaussian symplectic ensemble. Quarkonium 
system is time reversal and rotational invariant. There- 
fore in the case of chaotic motion its nearest-neighbor 
level spacing distribution is expected to be GOE-type. 
Now we carry out a detailed analysis of the spectra to 
establish the nearest-neighbor level spacing distribution 
function to determine the existence of quantum chaos in 
our system. 

The following values of the c and b quark masses and 
potential parameters are choosen in these claculations: 
For charmonium: 

m c = 1.486GeV , a s = 0.32, A 
0.2GeV 2 , OAGeV 2 , O.&GeV 2 . 
For bottomonium : 

m b = AMGeV , a s = 0.22, A 
OAGeV 2 , 0.6GeV 2 ,0.8GeV 2 . First 1000 eigenavleus 
are used to plott each histogram. In our claculations we 
vary only the principal quantum number n, keeping the 
orbital quantum number I constant. 

Figs. 1 and 2 show the level spacing distribution for 
charmonium and bottomium, respectively for Coulomb 
plus linear potential (k = 1) case, and in Fig.s 3 and 4 
the level spacing distributions are plotted for Coulomb 
plus harmonic oscillator potential (k = 2) case, at var- 
ious values of the confining potential parameter. It is 
clear that for the smaller values of the parameter A this 
distribution is close to GOE, that means the quantum 
system is becoming chaotic. It allows us to estimate the 
values of A where a transition from regular to chaotic 
behaviour will occur. Furthermore we find that for the 
same values of A the bottomonium motion is more reg- 
ular than that of the charmonium. Comparison of the 
level spacing distribution for cc and bb systems shows us 
that the transition from regular to chaotic motion for 
these systems occurs at different values of A. So, for cc 
system (Figs.l and 2) we have GOE-type distributions at 
A = 0.2GeV 2 and A = OAGeV 2 and Poissonian distribu- 
tion at A = 0.6GeV 2 , while the motion of bb system is still 
chaotic at A = 0.6GeV 2 (Figs.3 and 4). At A = 0.4GeV 2 
and A = 0.8GeV 2 the level spacing distribution for bot- 
tomonium is GOE and Poisson type, respectively. Thus 
it may be concluded from this comparison that the mo- 
tion of bb system is more chaotic than that of the cc 
system. 

Indeed, if one writes the Schrodinger equations for 
charmonium and bottomonium in the units m q = K = 1, 
the value of the confining potential for bottomonium is 
much smaller than that of charmonium. The chaotic na- 
ture of the motion in the quarkonium potential can be 
also understood from the following fact: formally the 
Hamiltonian for Coulomb plus harmonic oscillator po- 
tential is equivalent to the Hamiltonian of an atom in an 
uniform magnetic field, which is a nonintegrable system. 
It is well known, that the Hamiltonian of an atom in an 
uniform magentic field can be reduced to the Hamilto- 
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nian of coupled nonlinaer oscillators [10|, which exhibits 
chaotic dynamics. With this analogy the change in the 
confining potential parameters to some extent is equiva- 
lent to the change of the shape of billiard domain. 

IV. CONCLUSIONS 

Summarizing, in this work we have applied the tech- 
nique well-known in nuclear physics, atomic physics and 
other areas of physics to treat quantum chaos in heavy 
quarkonia, using for their description the Schrodinger 
equation with Coulomb plus confining potentials, a non- 
integrable system. It should be emphasized that this sys- 
tem is quite similar to the case of an atom in a constant 
magnetic field. It is well known that such a system is 
non-integralble and exhibits chaotic dynamics in both 
classical as well as quntum cases ^J. The analysis 
of the nearest level spacing distributions of these sys- 
tems shows, that for some smaller values of the confin- 
ing potential parameter the motion of the heavy quarko- 
nium is chaotic while for larger values of A the level spac- 
ing distribution is a Poisson type, that implies that the 
motion is becoming regular. Previously |24j the quan- 
tum chaos is related to the color screening parameter, 



while the present work treats the transition from regular 
to a chaotic motion related to a confining potential pa- 
rameters. The quarkonia can exhibit regular as well as 
chaotic dynamics depending on the value of the confining 
force. Study of the level statistics and quantum chaos in 
hadrons could play the same role in particle physics as 
it did in the case of nuclear physics. Therefore quantum 
chaology of hadrons could become one of the exciting top- 
ics in particle physics. More comprehensive analysis of 
the statistical properties of the heavy quarkonia spectra 
performed by solving the Dirac equation with Coulomb 
plus a confining potential is needed. This would establish 
the role played by relativistic effects in the approach to 
a level-spacing distribution. 



Acknowledgements 

The work of DUM is supported by NATO Science 
Fellowship of Natural Science and Engineering Research 
Council of Canada (NSERCC). The work of FCK is sup- 
ported by NSERCC. The work of URS is spupported by 
a Grant of the Uzbek Academy of Sciences (contract No 
33-02). 



[1] E.P. Wigner, Ann. Math. 53 36 (1951); ibid 55 7 (1952); 
ibid 62 548 (1955). 

[2] C.E.Porter, editor, Statistical Theories of Spectra (Aca- 
demic, New York. 1965); 

C.E.Porter and R.G.Thomas, Phys. Rev. 104 483 (1960). 
[3] M. L.Merita, Random Matrices and Statistical Theory of 

Energy Levels (Academic, New York. 1967); 

M.L.Mehta, Nucl.Phys. 18 395 (1960). 
[4] F.J.Dyson, J. Math. Phys. 3 157 (1962); ibid 3 166 (1962); 

ibid 3 1191 (1962). 
[5] T.A.Brody et al, Rev.Mod.Phys. 53 358 (1981). 
[6] B. Eckhardt, Phys. Rep. 163 207 (1988). 
[7] M.C.Gutzwiller, Chaos in classical and quantum systems. 

New York, Springer Verlag 1990. 
[8] O.Bohigas, M.J.Giannoni, and C.Schmit, Phys. Rev. Lett, 

52, 1 (1984). 

[9] S.H.Seligman, J.J.M.Verbaarschot, and M.R.Zirnbauer, 

J. Phys. A. 18 2751 (1985). 
[10] D.Delande and J.C.Gay, Phys.Rev.Lett., 57 2006 (1986). 
[11] H.Friedrich and D.Wintgen, Phys.Rep., 183 38 (1989). 
[12] T.Guhr, A.Mueller-Groeling, and H.A.Weidenmueller, 

Phys. Rep. 299 189 (1998). 
[13] ALEPH Collaboration, R.Barate, et al, Phys.Lett. B, 

425 215 (1998). 
[14] C.Quigg and J.L.Rosner , Phys.Rep., 56 167 (1979). 
[15] W.Lucha, F.F.Schoberl and D.Gromes, Phys.Rep., 200 

128 (1991). 



[16] X.Q.Zhu, F.C.Khanna, R.Gourishankar and R.Teshima, 

Phys.Rev.D, 47 1155 (1993). 
[17] S.N.Mukherjee et al., Phys.Rep. 231 203 (1993). 
[18] M.A.Halasz and J.J.M.Verbaarschot, Phys.Rev. Lett., 74 

3920 (1995). 

[19] . A. Berg, H. Markum, and R. Pullirsch Phys. Rev. D 59 

097504 (1999) 

[20] Vladimir Pascalusta, hep-ph/0201040 
[21] B.A.Berg, et.al hep-lat/0007008 

[23] E.Bittner, H.Markum and R.Pu llirsch, \hep-lat/0110222\ 

[23] D. Bukta, G. Karl, B. G. Nickel chao-dyn/9910026 

[24] Jian-zhong Gu et al, Phys.Rev. C, 60 035211 (1999). 

[25] L.A.Bunimovich, Commun.Math.Phys., 65 259 (1979). 

[26] M.V.Berry, Eur.J.Phys., 2 91 (1981) 

[27] G.M.Zaslavsky, Stochastic Behaviour of Dynamical Sys- 
tems. New York, Harwood. 1985. 

[28] Francisco M. Fernandez and Eduardo A.Castro, J.Chem. 
Phys., 79 321 (1983). 

[29] D.U.Matrasulov, M.M.Musakhanov and T.M.Morii, 
Phys.Rev. C, 61 045204 (2000) 

[30] D.U.Matrasulov, F.C.Khanna and H.Yusupov J. Phys. G. 
29 473 (2003) 

[31] L.D. Landau and E.M.Lifshitz, Quantum Mechanics, Ox- 
ford, 1965. 

[32] M.V.Berry Proc. Roy.Soc. London. 400 229 (1985) 



This figure "figl.gif" is available in "gif" format from: 



http://arXiv.org/ps/hep-ph/0306214vl 



This figure "fig2.gif" is available in "gif" format from: 



http://arXiv.org/ps/hep-ph/0306214vl 



This figure "fig3.gif" is available in "gif" format from: 



http://arXiv.org/ps/hep-ph/0306214vl 



This figure "fig4.gif" is available in "gif" format from: 



http://arXiv.org/ps/hep-ph/0306214vl 



